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Abstract. Recent analysis of the helioseismic observa- 
tions indicate that the previously observed surface tor- 
sional oscillations extend significantly downwards into the 
solar convection zone. 

In an attempt to understand these oscillations, we 
study the nonlinear coupling between the magnetic held 
and the solar differential rotation in the context of a mean 
held dynamo model, in which the nonlincarity is due to the 
action of the azimuthal component of the Lorentz force of 
the dynamo generated magnetic field on the solar angu- 
lar velocity. The underlying zero order angular velocity is 
chosen to be consistent with the most recent helioseismic 
data. 

The model produces butterfly diagrams which are in 
qualitative agreement with the observations. It displays 
torsional oscillations that penetrate into the convection 
zone, and which with time migrate towards the equator. 
The period of these oscillations is found to be half that 
of the period of the global magnetic fields. This is com- 
patible with the observed period of the surface torsional 
oscillations. Inside the convection zone, this is a testable 
prediction that is not ruled out by the observations so far 
available. 
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1. Introduction 

An important feature of the solar convection zone is the 
presence of differential rotation in the form of a decrease 
in angular velocity from equator to the pole. This has been 
observed both in the surface layers (e.g. Snodgrass 1984) 
and deeper in the convection zone, as inferred from helio- 
seismic measurements (e.g. Thompson et al. 1996). Fur- 
thermore, the differential rotation on the surface has been 
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observed to vary with time (e.g. Howard & LaBonte 1980; 
Snodgrass, Howard, & Webster 1985). These so called tor- 
sional oscillations, which have periods of about 11 years, 
manifest themselves in the form of four alternating lat- 
itudinal bands of slightly faster and slower than average 
zonal flows which migrate towards the equator in about 22 
years. These oscillations have also been confirmed by the 
analysis of the helioseismic data from the Solar and He- 
liospheric Observatory (SOHO) spacecraft for the present 
solar cycle (Kosovichev & Schou 1997; Schou et al. 1998). 

Recent analysis of the helioseismic data, from both the 
Michelson Doppler Imager (MDI) instrument on board 
the SOHO spacecraft and the Global Oscillation Network 
Group (GONG) project, has also produced evidence that 
this banding signature is not merely a surface feature, but 
extends into the convection zone, to a depth of at least 
8 percent in radius (Howe et al. 2000). These authors 
present data on departures of the reconstructed rotation 
rate from its temporal averages - the residuals - as a func- 
tion of latitude at several target depths, which behave in 
a manner similar to the migration of sunspots during the 
solar cycle (the 'butterfly diagram'). This finding is also 
supported by the analysis of Antia & Basu (2000), who 
use different data sets from GONG and independent in- 
version techniques. The time-base of these observations is 
only a few years, less than a complete solar cycle. 

These torsional oscillations are thought to be produced 
as a consequence of the nonlinear interactions between 
the magnetic fields and the solar differential rotation. A 
zero order 'mean' differential rotation is assumed to be 
maintained by the Reynolds stresses of the turbulence; 
this can be included as a constant 'background' angular 
velocity, or explicitly parametrized, e.g. by the so called 
A-effect (e.g. Rudiger 1989). Attempts have been made 
to explain the surface oscillations in terms of the effects 
of the Lorentz force exerted by the large scale magnetic 
field on the azimuthal velocity field (e.g. Brandenburg & 
Tuominen 1988), or as a consequence of the 'quenching' of 
the turbulence-dependent quantities by the magnetic field 
(Kitchatinov et al. 1999). 

Here we study this nonlinear coupling in the context 
of a two dimensional axisymmetric mean field dynamo 
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model, in a spherical shell, in which the only nonlinearity 
is the action of the azimuthal component of the Lorentz 
force of the dynamo generated magnetic field on the solar 
angular velocity. Obtaining torsional oscillations in this 
way is also of interest in view of the fact that the Lorentz 
force is of second order in the magnetic field, thus natu- 
rally leading to the excitation of hydrodynamical oscilla- 
tions of about half the period of the magnetic oscillations. 

In the next section we introduce our model. Section 3 
contains our results and section 4 gives a brief discussion. 

2. The model 

Here we shall assume that the gross features of the large 
scale solar magnetic field can be described by a mean field 
dynamo model, with the standard equation 



$B 

— = V x (u x B + aB - 77 V x B) 



(1) 



where u = vcf> — 5V77, the term proportional to V77 rep- 
resents the effects of turbulent diamagnetism, and where 
the velocity field is taken to be of the form 



v = vo + v 



(2) 



where vq = f2o?"sin#, Qq is a prescribed underlying rota- 
tion law and the component v' satisfies 

dv' _ (V x B) x B • 
~dt 
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where D 2 is the operator ■§ pl + f-§p + ^jr^ 
-p-s) and fiQ is the induction constant. The assumption of 
axisymmetry allows the field B to be split into toroidal 
and poloidal parts, B = Bp + Bp = B(f> + V x A<j>, and 
results in Eq. ([!]) yielding two scalar equations for A and 
B. Nondimcnsionalizing in terms of the solar radius Rq 
and time Rq/tjo, where 770 is the maximum value of 77 
and letting = fTfi, a = a^a, 77 = 77077, B = t3oB 
and v' = Q*Rqv', results in a system of equations for 
A, B and v' , with the dynamo parameters R a = a^RQ/rjo, 
R u = Q*Rq/ti 0i P r = Vo/r] , and 77 = 77/770, where tt* is 
the solar surface equatorial angular velocity (see Moss & 
Brooke 2000 for details) . Here vq and 770 are the turbulent 
magnetic diffusivity and viscosity respectively and P r is 
the turbulent Prandtl number. The density p is assumed 
to be uniform, and stress free boundary conditions ensure 
angular momentum conservation. 

These equations were solved using the code and bound- 
ary conditions described in Moss & Brooke (2000), over 
the range tq < r < 1, < 9 < 7r, with uniform spac- 
ing in both r and 9. The computational domain is the 
region r$ = 0.64 < r < 1; with the solar convection zone 
proper being thought to occupy the region r > 0.7, the 
region r < r < 0.7 can be thought of as an overshoot 
region/tachocline. In the following simulations we used 
a mesh resolution of 61 points uniformly distributed ra- 
dially and 101 points uniformly distributed latitudinally 



(over < 9 < tt) , but test runs were carried out at higher 
spatial resolutions. 

In this investigation, we took £!o to be given in 0.64 < 
r < 1 by an interpolation on the MDI data obtained from 
1996 to 1999 (Howe et al. 2000), depicted in Fig. 1. For 
a we took a = a r (r)f(9), where f(9) = sin 2 6 l cos0 (cf. 
Rudiger & Brandenburg 1995) and 
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0.7 < r < 0. 



(4) 



with cubic interpolation to zero at r = ro and r = 1, with 
the convention that a r > and R a < 0. Also, in order 
to take into account the likely decrease in the turbulent 
diffusion coefficient 77 in the overshoot region, we allowed 
a simple linear decrease from f] = 1 at r = 0.8 to f) = 0.5 
in r < 0.7. 

We monitor the time evolution of the total magnetic 
energy E and the global parity of the magnetic field, de- 
fined as P = gsznpt i where 5 and A refer respectively to 
the parts of the magnetic field that have symmetry or an- 
tisymmetry with respect to the equatorial plane (see also 
Brandenburg et al. 1989). Thus P = +1 and —1 corre- 
spond to symmetric and antisymmetric fields respectively. 

3. Results 

We calibrated our model so that near marginal excita- 
tion the cycle period was about 22 years. This determined 
Ru = 6 x 10 4 , corresponding to 770 ~ 2.5 x 10 11 cm 2 s _1 , 
given the known values of Q* and Rq . The first solutions 
to be excited in the linear theory are odd parity (P = —1) 
limit cycles, which in this case have marginal dynamo 
number R a ~ —3.12. The even parity (P = +1) solutions 
are also excited at similar marginal dynamo numbers of 
R a ~ —3.16. We considered two values of the Prandtl 
numbers, P r = 0.1 and P r = 1.0. 

With these parameter values, we found that this 
model, with underlying zero order angular velocity cho- 
sen to be consistent with the recent (MDI) helioseismic 
data (Fig. 0), is capable of producing butterfly diagrams 
which are in qualitative agreement with the observations. 
An example of this is depicted in Fig. ||. The polar feature 
is rather too strong - we have checked that this can be 
rectified by modifying slightly the spatial dependence of 
a, by for example choosing f(9) = sin 9cos9. 

We also found that this model successfully produced 
torsional oscillations in the convection zone, similar to 
those deduced from recent helioseismic data. To compare 
our model with these results, we have plotted in Figs. || 
and ^ the variations of the rotation rate with latitude and 
time ('butterfly diagrams') to reveal the migrating banded 
zonal flows. In these models, the basic magnetic field par- 
ity is odd (P = -1). 

For the sake of comparison with observational results, 
we have also plotted in Fig. [| the evolution of the residual 
rotation rate with time, at radius 0.84 and latitude 60 
degrees. 
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Fig. 1. Isolines of the time average of the angular velocity 
of the solar rotation, obtained by inverse techniques using 
the MDI data (Howe et al. 2000). Contours are labelled in 
units of nHz. 




Time (years) 

Fig. 2. Butterfly diagram of the toroidal component of 
the magnetic field B at R = 0.95-Rq. Dark and light 
shadows correspond to positive and negative respec- 
tively. Parameters values are R a = —3.2, P r = 1.0 and 
R u = 6 x 10 4 . 

As can be seen, consistent with the observations, in 
each hemisphere there are alternating latitudinal bands, 
with the width of approximately 10 degrees, of slightly 
faster and slower than average zonal flows. These migrate 
towards the equator in about 22 years, and extend deep 
into the convection zone. The amplitudes of these oscilla- 
tions increase with depth below the surface and depend 
on the parameters of our model, in particular the Prandtl 
number. For P r = 1.0, these amplitudes range from about 
0.07 nHz at the surface to more than 0.4 nHz towards the 
bottom of the convection zone, somewhat lower than, but 
in principle compatible with, the results of Howe et al. 
2000. The torsional oscillations present in our model have 
periods half that of the period of the global magnetic field 




Fig. 3. Variation of rotation rate with latitude and time 
from which a temporal average has been subtracted to re- 
veal the migrating banded zonal flows, at R = O.95i? 
(top) and R = 0.84i? Q (bottom). Darker and lighter re- 
gions represent positive and negative deviations from the 
time averaged background rotation rate. Parameters val- 
ues are R a = -3.2, P r = 0.1 and R u = 6 x 10 4 . 

which is compatible with the observed period of the os- 
cillations at the surface and consistent with the observed 
behaviour inside the convection zone. 

The strictly odd parity models presented have an equa- 
torial feature in the 'butterfly diagrams' for the velocity 
perturbations which does not appear to be present in the 
current inversions of the observational data. This becomes 
weaker as Prandtl number increases. We also note in pass- 
ing that these figures are almost identical for cases where 
even parity solutions are found, except that this feature is 
then absent. We further point out that the large scale solar 
magnetic field is probably of mixed global parity (predom- 
inantly odd) - sec, e.g., Pulkinncn et al. (1999). 

4. Discussion 

We have studied a solar dynamo model, calibrated to have 
the correct cycle period, with a mean rotation law given 
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Fig. 4. Variation of rotation rate with latitude and time 
from which a temporal average has been subtracted to 
reveal the migrating banded zonal flows, at R — O.95i? 
(top) and R = 0.84i? Q (bottom). Parameters values are 
R a = -3.2, P r = 1.0 and i? w = 6 x 10 4 . 



We have shown that a nonlinear coupling between the 
magnetic field and the solar differential rotation, where 
the nonlinearity is due to the action of the azimuthal com- 
ponent of the Lorentz force, is capable of producing tor- 
sional oscillations, with period of about 11 years, which 
penetrate into the convection zone and which migrate to- 
wards the equator in about 22 years. The period of these 
oscillations is about half that of the period of the global 
magnetic fields. This is in agreement with the observed 
period of the torsional oscillations at the surface. For the 
oscillations inside the convection zone, this is a testable 
prediction, not contradicted by the current helioseismic 
observations, which so far extend over an interval less than 
a solar cycle. We note also that solutions with even par- 
ity (P = +1) show slightly larger amplitudes without the 
equatorial feature. 

The current inversions of the helioseismological data 
seem to suggest that the torsional oscillations largely dis- 
appear below about R — O.9i?0, in contrast to our model 
oscillations. However, there are uncertainties in these in- 
versions, specially at the deeper levels. At the same time 
our dynamical model is oversimplified and substantial im- 
provements can be made. In particular, the predicted am- 
plitudes for torsional oscillations in our model are likely to 
be affected by our assumption of uniform density in Eq. 

. Nevertheless, we find it interesting that some of the 
major features in the torsional oscillations can be readily 
reproduced. A more detailed study of these oscillations, 
including more realistic density profiles, is in progress. 
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Fig. 5. Variation of rotation rate at radius 0.84i? Q and 
latitude 60 degrees. Also depicted is the (nondimensional- 
ized) toroidal field at the same radius and latitude (dashed 
line), which shows double the period. Parameters values 
are as in Fig. 4. 



by recent helioseismic observations. This model produces 
butterfly diagrams in qualitative agreements with those 
observed. 
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